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In the context of the dynamical mean field theory of the Hubbard model, we study the behavior of 
the compressibility near the density driven Mott transition at finite temperatures. We demonstrate 
this divergence using dynamical mean field theory and Quantum Monte Carlo simulations in the 
one band and the two band Hubbard model. We supplement this result with considerations based 
on the Landau theory framework, and discuss the relevance of our results to the a-y endpoint in 
Cerium. 



PACS Numbers: 71.30. +h, 71.10.Fd, 71.27.+a 

The Mott transition, namely the metal-insulator tran- 
sition (MIT) driven by electron-electron interactions [Q, 
is a fascinating phenomenon realized experimentally in 
many compounds such as V2O3 and Ni(Se,S)2 |^]. The 
Mott transition concept is also relevant to elements in 
the lanthanide and actinide series Viewed from a 
broader perspective, the Mott transition problem forces 
us to develop tools to describe materials where the elec- 
tron is not fully described by either a real space picture 
or a momentum space picture, and continues to spur ad- 
vances in many body and electronic structure methods. 

On the theory side, the Hubbard model is the sim- 
plest Hamiltonian that captures some of the essential 
physics of the transition. It has been intensively stud- 
ied in one dimension, but in this limit no finite temper- 
ature phase transitions can take place. In recent years, 
theoretical progress has been made in the understand- 
ing of the Mott-Hubbard transition using the dynamical 
mean-field theory (DMFT) g. In this framework, Mott 
transition points can be viewed as bifurcation points of 
a functional of the local Greens function G, or of its 
associated variable, the Weiss field which describes the 
local environment of a correlated site. 

The case of the correlation strength (U) driven MIT 
at half filling, is now well understood. At temperature 
T = there are two bifurcation points, one denoted by 
Uci{T — 0) where the insulating solution disappears, and 
the other denoted by Uc2{T = 0) where the metallic state 
disappears in a fashion reminiscent to the Brinkman Rice 
scenario [0 . It was found that in the U-T phase diagram 
of the frustrated Hubbard model there is a first-order 
MIT line that ends in a finite temperature second- 
order critical point {Tmit, Umit) which has the charac- 
ter of a regular Ising bifurcation with a rapid variation 
of the susceptibility connected to the double occupancy 
0,|l|. At higher temperatures the Uc2iT) and Uci{T) 
lines become crossover lines, which have a well defined 
experimental significance discussed in |p^ . 

The zero temperature aspects of the doping driven 
MIT were studied in ||l^ . It was shown that there are two 
solutions in an area bound by the curves jJci, where the 



insulating solution disappears, and fj,c2, where the metal- 
lic state disappears. The finite temperature aspects of 
the doping driven Mott transition have not been investi- 
gated so far. This is the subject of this paper. We will 
not consider the effects of long range order. 

Our main interest is the behavior of the charge com- 
pressibility near the doping driven Mott transition in 
the paramagnetic phase at finite temperatures. Fu- 
rukawa and Imada pointed out that the compress- 
ibility diverges at the density driven Mott transition in 
2-dimensions at T = 0. This behavior has been also 
observed on other models of correlated electron systems 
such as the t-J model indicating that this phenomenon is 
quite general [|l5|. Simple models of the Mott transition, 
such as the Gutzwiller approximation or the slave-boson 
approach predict a finite compressibility psf . It is very 
important to understand the physical origin of this re- 
sult, and to see if it is realized in the DMFT solution of 
the Hubbard model. The previously investigated bifurca- 
tion points within the DMFT, have either a finite charge 
compressibility, such as in the T = density driven Mott 
transition, or a vanishing charge compressibility, as in the 
T — correlation driven transition. 

Our study of the neighborhood of the Mott transition 
endpoint is relevant to materials which have a finite tem- 
perature isostructural phase transition such as Cerium. 
This can be seen by generalizing the derivation of the 
Landau free energy of [|ll[. Near the transition, the free 
energy of more complicated models would have the same 
form as that studied in our paper. 

We start with the m-band degenerate Hubbard model: 
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The first term describes the hopping between nearest 
neighbors {ij) on a lattice with coordination number z. 
m,m' — 1,2 are the band indices and a =|, I labels the 
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spin index. The parameter U is the energy cost associ- 
ated with having a double occupancy on each site, t and 
U are assumed to be independent of the band indices and 
we concentrate on the paramagnetic phase. In the Umit 
of infinite dimensions, z — > oo, this model can be mapped 
onto a single-impurity Anderson model (SIAM) supple- 
mented by a self-consistency condition Q. The DMFT 
equation of the model reads 



(2) 



where uj„ are the fermionic Matsubara frequencies, Gma 
is the impurity Green's function and A is the hybridiza- 
tion function of the SIAM. Here, a denotes the control 
parameter such as temperature T, Coulomb repulsion U, 
chemical potential /i, etc. We adopt a semicircular den- 
sity of states, /9(e) = — {e/DY where the half- 
bandwidth D = 2t — 1 la our unit of energy. 

We now propose a schematic phase diagram (Fig. ^ for 
the degenerate Hubbard model. The figure shows cross 
sections of the T-^, plane for various values of the interac- 
tion U . In the particle-hole symmetric case in the 1-band 
model the peaks have equal height. The shaded regions 
indicate coexistence between metallic and insulating so- 
lutions. The presence of a coexistence region and a first 
order phase transition in 2-band model was apparent in 
earlier Monte Carlo calculations This also follows 

from the general form of the Landau functional of this 
problem Q , which has the same form for the 1-band and 
and the multiband situations. The finite temperature 
and general chemical potential aspects of this problem 
had not been discussed before in the literature. 




FIG. 1. Schematic phase diagram for the degenerate Hub- 
bard model. The cross sections shown are on the T-p plane 
for different values of U . jici (the heavy line) and Ud are 
the chemical potential and interaction respectively at which 
the insulating solution gets destroyed. /Xc2 (the dotted line) 
and Uc2 are the chemical potential and interaction at which 
the metallic solution gets destroyed. Umit is the value of the 
interaction at which the MIT endpoint takes place. 



On solving the DMFT equation (^ iteratively using 
quantum Monte Carlo (QMC) methods, we find that at 
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FIG. 2. Particle occupation 5 — (n)-l/2 per spin as a func- 
tion of the shifted chemical potential /i — (7/2 for different 
temperatures T. The lower panel shows QMC data obtained 
at r = 1/40, 1/30, 1/20, 1/16 (right to left) which are above 
Tmit- The upper panel shows similar data at temperatures 
T — 1/64, 1/50 which are below Tmit- The left arrow indi- 
cates that starting at high doping, a metallic solution can be 
followed down to half-filling. Also, starting from an insulat- 
ing state at half-filling, one can follow the coexisting solution 
(right arrow) by increasing /i up to a sudden jump in S. 



integer filling, for smaller values of t/ (< Uci in Fig. |^) 
the solutions show metallic behavior while for large val- 
ues of J7 (> Uc2 in Fig. |l|) the solutions show insulating 
behavior. In Fig.|| we show the results for the doping 
(per spin) S = (n) — 1/2 as a function of the chemical po- 
tential /i obtained from QMC calculations with At = 0.5 
in the single band case. The interaction value U = 2.46 
is inside the small region of coexisting solutions when 
the system is at half- filling |T^. In the lower panel we 
show the results at several temperatures above the criti- 
cal temperature Tmit- As the temperature is decreased 
the curves develop a sigmoidal shape, which is a hall- 
mark of the approach to a second order critical point in 
Landau theory, as in the familiar Ising mean-field model. 

The results for the total occupation number n as a 
function of /i in the 2-band model are shown in Figj^. 
The QMC calculation was carried out at At = 0.25 and 
U = 3.0. At higher temperatures, only one solution is 
present, whereas on lowering the temperature to T = 
1/40 both metallic and insulating solutions coexist. 

To obtain coexisting solutions below Tmit we first 
start at high doping in a metallic state as indicated by the 
left arrow in the upper panel of Fig.|[ For this, we use a 
seed from low U (< Uci in Fig.|^) and continuously evolve 
towards integer filling through the upper branch solution 
which always remains metallic with finite compressibil- 
ity. We use a self-consistent solution as new seed for the 
successive calculations at smaller /i. On the other hand. 
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FIG. 3. Particle occupation n function of [i, for different 
temperatures, in the 2-band model at U=3.0. In the upper 
panel, the temperatures are above the critical temperature 
while in the lower panel, we see coexistence at T < Tmit- 



a second solution is also present at half-filling in Fig.g 
(or near n = 1 in Fig.^) corresponding to an insulting- 
like state, obtained by using a seed from large U (> Uc2 
in Fig. 1^). This state is essentially incompressible as n 
almost remains constant while increasing /x. This can be 
continuously followed as fj, is increased until the even- 
tual jump of n towards the unique solution present at 
the higher values of ^. This procedure successfully de- 
termined the location of the coexisting region ||l8| . 

Locating these parameters in the schematic phase dia- 
gram, we find that the data for the 1-band model below 
Tmit in Fig. || corresponds to the region near the right 
peak at Umit < U < Uc2 in Fig. |l]. As we move towards 
(5 = in the metallic phase, we cross the coexistence re- 
gion and reach the ^c2 line at T = 1/50. Upon lowering 
T to 1/64 we can approach 5 — Q without crossing the /Zc2 
line. In the 2-band model, it was found that at T = 1/40 
the insulating solution disappears as we move away from 
n — 1 indicating the the ^ci line has been crossed indi- 
cating that the data in Fig. || lies near the right peak at 
Umit <U <Uc2 in Fig. g. 

The region with two coexisting solutions has two dif- 
ferent values of n for given T and /i. These two coexisting 
solutions have different free energies and the actual ther- 
modynamic state of the system will be that of minimum 
energy. Therefore, a jump in the particle occupation is 
predicted at a first-order line. The actual determination 
of this line, thus implies a precise calculation of the free 
energy which lies outside the scope of the present work. 

From the n vs. fi curves in Figs. I and I we com- 
puted the numerical derivative of the particle number 



FIG. 4. The inverse compressibility k,~^ , obtained by the 
numerical differentiation of n with respect to fi, at a constant 
doping as a function of T. Top: 1-band model at [/ = 2.46 
with doping 5 = 0.002. Bottom: 2-band model at U = 3.0 at 
n = 1.003. The solid lines are a guide to the eye. 



with respect to the chemical potential to obtain charge 
compressibility k. The results for function of the 

temperature above Tmit (Fig.^ indicate that k^^ ^ 
as we pass through the MIT. 

We also like to mention that in our simulations we ob- 
served characteristic effects of enhanced fluctuations and 
critical slowing down as the MIT is approached. Hence, 
simulations have to be done with extreme care, appro- 
priately choosing the seeds for the iterative process and 
gathering many sample points for accurate statistics. In 
practice we use ~ 10^ MC sweeps and run our calcula- 
tions in an 8 node parallel cluster. As we approach the 
critical points, we used up to a few hundred iterations to 
obtain converged solutions. 

We now complement our numerical results with Lan- 
dau theory arguments. The mean-field equation (H) can 
be obtained by differentiating the Landau functional 



i^imp[A] 



(3) 



with respect to the hybridization A(iti;„) of the SIAM, 
which has the meaning of a Weiss field. i^imp[A] is the 
free energy of the SIAM in presence of a hybridization. 
The Green's function of the SIAM is G(zw„)[A,a] = 
(l/2T)(5i^imp/(5A. This Landau approach was used to 
describe the Mott transition at half filling [|lO|,0 . 

As discussed in Ref. |ll|, the finite temperature Mott 
transition is a regular bifurcation point. On differentiat- 
ing the Landau functional twice we get a matrix of the 
form —S„r 
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acquires a zero mode. 

We now make a small change in chemical potential 
around the critical point and expand the mean field 
equation(|^) to first order in Sa = (/i — /ic), SA — 
A{ac + 6a) - A{ac)- We get, 
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From the lattice occupation at any site which is 
identical to the local impurity occupation is related to 
the hybridization by (n) = {2T/t'^)J^^A{iujn)- Thus, 
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Clearly, unless the derivative in the rhs of (g) is exactly 
orthogonal to the zero mode of the matrix M, the bifurca- 
tion condition leads to the singular behavior of the com- 
pressibility. The QMC studies shown above demonstrate 
that this orthogonality does not occur. The divergence of 
the compressibility strengthens the liquid-gas picture of 
the Mott transition presented by Castellani et. al. 

From the experimental viewpoint, we believe that our 
results highlight important aspects of the physics of the 
a-"f transition in Cerium. The detailed description of 
the non-universal aspects of this material requires more 
quantitative studies or more elaborate models such as 
those studied by Held et. al. [|2^. However, the func- 
tional describing these complicated models near the finite 
temperature Mott transition would reduce to the one un- 
derlying the equations we studied. Hence it is plausible 
that the behavior of the compressibility that we identified 
in our basic model applies to the a-7 transition of Cerium 
as well. The divergence of compressibility in the Cerium 
a-j transition therefore has an electronic origin and can 
be understood from model calculations without involv- 
ing lattice degrees of freedom (which would renormalize 
values of the critical points without changing qualitative 
features) . In our data, we find a decrease in compressibil- 
ity during the transition from the insulating to the metal- 
lic phase which is similar to what has been measured by 
Beecroft and Swenson in ||2l[. According to the results 
in ]2l|] , the compressibility in the low-pressure 7 phase 
(corresponding to the insulating region) is larger than in 
the high pressure a phase (corresponding to the metallic- 
like state) at low temperatures, and this effect becomes 
less noticeable at higher temperatures. Our QMC stud- 
ies have shown that there is a region of coexistence where 
the metallic phase is extended into the stability region of 
the insulating phase, corresponding to the hysteresis ob- 
served in the data in |21|. This was also found in X-ray 
diffraction studies on Cerium [ p2[ , where the 7 phase was 
quenched into into the region of stability of the a phase 
near the phase boundary. 

In summary, we presented a careful QMC numerical 
solution of the doping driven Mott transition in the limit 



of large lattice connectivity. Our study unveils that the 
divergence of the compressibility is a generic feature of 
the confluence of the flnite temperature Mott endpoint. 
We understood this divergence in terms of a simple ar- 
gument based on Landau theory, which indicates that 
these results are more general than the specific models 
for which the numerical studies were carried out. Our 
results were found to be relevant to the physics of the 
Cerium a-7 transition. 
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